A cell functional minimization scheme for parabolic problem

We are interested in a robust and accurate finite volume scheme for 2-D parabolic problems derived from the cell functional minimization approach. The scheme has a local stencil, is locally conservative, treats discontinuity rigorously and leads to a symmetric positive definite linear system. Since the scheme has both cell centered unknowns and cell edge unknowns, the computational cost is an issue and a parallel algorithm is then suggested based on nonoverlapping domain decomposition approach. The interface condition is of the Dirichlet–Robin type and has a parameter λ. By choosing this parameter properly, the convergence of the iteration process could be sped up. Numerical results for linear and nonlinear problems demonstrate the good performance of the cell functional minimization scheme and its parallel version on distorted meshes.     

Induced polarization forward modelling using finite element method and the fractal model

The modeling of induced polarization data from 2D and 3D geological medium was developed. We applied the finite element method and the fractal model to complex resistivity of rocks as an intrinsic electrical property of geological medium. The LU factorization and pre-conditioned complex bi-conjugated gradient method were applied to solve the complex linear equation system. The modeling implemented in this work was efficient and prompt. In addition, the induced polarization survey using the fractal model as intrinsic property of the medium is a promising method for environmental investigation.     

Velocity–vorticity–helicity formulation and a solver for the Navier–Stokes equations

For the three-dimensional incompressible Navier–Stokes equations, we present a formulation featuring velocity, vorticity and helical density as independent variables. We find the helical density can be observed as a Lagrange multiplier corresponding to the divergence-free constraint on the vorticity variable, similar to the pressure in the case of the incompressibility condition for velocity. As one possible practical application of this new formulation, we consider a time-splitting numerical scheme based on an alternating procedure between vorticity–helical density and velocity–Bernoulli pressure systems of equations. Results of numerical experiments include a comparison with some well-known schemes based on pressure–velocity formulation and illustrate the competitiveness on the new scheme as well as the soundness of the new formulation.     

半线性抛物方程的时空有限元方法的误差估计

<正>1引言本文考虑如下半线性抛物方程(?)其中Ω∈R~2.函数f(u):C→C满足:(1)|f(u)|≤c|u|(?)u∈C(Ω)(2)Lipschitz条件,即     

Product Inhibited Enzymatic Reactions on the Rotating Disk Electrodes

A model of cyclic staircase voltammetry of enzymatic reactions on the rotating disk electrode is developed and the inhibition by the product is investigated. The responses are mostly time – dependent current – potential curves with maxima in the forward but not in the backward branch. If the maxima appear in both branches, the inhibition is potential – dependent either directly or indirectly. The source of the second maximum is explained.     

A Petrov-Galerkin method for equal width equation

The equal width equation is solved numerically by Petrov-Galerkin method using linear hat function and quadratic B-spline function as trial and test functions respectively. Product approximation has been used in this method. A linear stability analysis of the scheme is shown to be conditionally stable. Test problems including the single soliton and the interaction of solitons are used to validate the suggested method, which is found to be accurate and efficient. The Maxwellian initial condition pulse and the development of an undular bore are also studied.     

基于后向散射模型的蒸汽湿度测量及仿真

为了更加准确地测量出汽轮机末级蒸汽湿度,为锅炉给水量提供依据,促进汽轮机内部除湿技术的发展,采用了后向角散射法的思路设计了湿蒸汽参数测量模型,并对模型所需物理量进行了详细的计算,为后向散射法测量蒸汽湿度的研究提供参考。通过对模型进行仿真模拟以及实验研究,得到了后向散射光光强数据。通过对比仿真模拟数据与实验数据,验证了该模型的可行性与有效性,对光散射法测量蒸汽湿度的实验分析以及工程实践具有重要意义。     

Wavelength selection of bidirectional laser transmission based on Monte Carlo method

The laser detection technology in uncertain and dynamic environments is of utmost importance in many fields. A model of transient radiative transfer of bidirectional path laser based on Monte Carlo method is developed to investigate the optimum wavelength of active detector at complex atmospheric conditions. The radiative parameters of atmosphere are calculated by HITRAN database and Mie theory at several typical atmospheric conditions including the standard atmosphere, urban aerosol, and radiation fog. Transmission characteristics for five spectral bands at the above atmospheric conditions are calculated by this model. The optimal transmission ability occurred in bands 0.2–0.5, 1.4–1.6, and 0.75–1.25 μm on the condition of standard atmosphere, urban aerosol, and radiation fog, respectively. All results provide effective reference and basic support for choosing the optimal spectral band for active detection.     

Finite size effects in stimulated laser pair production

We consider stimulated pair production employing strong-field QED in a high-intensity laser background. In an infinite plane wave, we show that light-cone quasi-momentum can only be transferred to the created pair as a multiple of the laser frequency, i.e. by a higher harmonic. This translates into discrete resonance conditions providing the support of the pair creation probability which becomes a delta-comb. These findings corroborate the usual interpretation of multi-photon production of pairs with an effective mass. In a pulse, the momentum transfer is continuous, leading to broadening of the resonances and sub-threshold behaviour. The peaks remain visible as long as the number of cycles per pulse exceeds unity. The resonance patterns in pulses are analogous to those of a diffraction process based on interference of the produced pairs. We finally comment on the dependence of the peak positions, and in turn the effective mass, on the pulse shape.     

Modified characteristic finite difference fractional step method for moving boundary value problem of percolation coupled system

For the coupled system with moving boundary values of multilayer dynamicsof fluids in porous media,a characteristic finite difference fractional step scheme appli-cable to the parallel arithmetic is put forward.Some techniques,such as the change ofregions,the calculus of variations,the piecewise threefold quadratic interpolation,themultiplicative commutation rule of difference operators,the decomposition of high orderdifference operators,and the prior estimates,are adopted.The optimal order estimatesin the l2norm are derived to determine the error in the approximate solution.This nu-merical method has been successfully used to simulate the flow of migration-accumulationof the multilayer percolation coupled system.Some numerical results are well illustratedin this paper.